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CLAIMS 

Claim 1 (original): Method for estimating the seismic illumination fold \(x,p) 
in the migrated 3D domain at least one image point x , for at least one dip of vector p , 
wherein the illumination fold 1(1, p\sj) for each (source J, receiver 7) pair in the 
seismic survey is estimated, by applying the following steps: 

- determination of the reflection travel time t r (x r (/?);?,r) from the source s 
to the specular reflection point x r on the plane reflector passing through the image point 
x and perpendicular to the dip vector p and then returning to the reflector 7 ; 

starting from the diffraction travel time t d {x\sj) from the source J to the 

said image point x and then returning to the reflector 7 ; 

- incrementing the said illumination fold l(x,p;s,7) related to the said (source 
J , receiver 7 ) pair as a function of the difference between the diffraction travel time 
t d {x\ s , 7) and the reflection travel time t r (x r (/?); s , 7 ) . 

Claim 2 (original): Method according to claim 1, comprising the step of 
summating each of the said illumination folds I( J, /?; s, 7) related to a (source s , receiver 

7 ) pair so as to determine the total illumination fold I( J, p) = ^ I(3c, p; 5", 7 ) . 

Claim 3 (previously amended): Method according to claim 1, wherein, during the 
incrementing step, the illumination fold I(x, p>s,7) is incremented using an increment 

function \{t d j r \sj) according to 1(1, p) = I(x,p) + i(^/ r ;5,r), the said increment 

function taking account of the difference between the diffraction travel time t d (x\sj) 

and the reflection travel time t r (x r {p)\ s, 7) . 

Claim 4 (original): Method according to claim 3, wherein the increment 
function i is a function of the seismic wavelet s(t). 
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Claim 5 (original): Method according to claim 4, wherein the increment 
function i is expressed as a function of the derivative of the seismic wavelet s(t) 
according to: 

i (t d >t r ;s, r)= s( t d (x; 7) - t r (x r (p)\ sj). 

Claim 6 (original): Method according to claim 4, wherein the increment 
function i is expressed as a function of the derivative s (t) of the seismic wavelet s(t) with 
respect to time according to: 
i (t d ,/,;?, 7) =( t d (x; 7) - t r (x r {p)\ 7) . 

Claim 7 (original): Method according to any one of claims 3 to 6, in which an 
a priori correction w (J, s, r) of the illumination fold is taken into account by migration, 
comprising the step of incrementing the illumination fold 1(1, p\sj) related to a (source 
s , receiver r ) pair by i (t d t r \sj). w( J; 5, 7) . 

Claim 8 (previously presented): Method according to claim 1, wherein the 
determination step includes the second order Taylor series development of the diffraction 
travel time t d (x; 5, r ) around the image point x : 

t d (x; r ) = t d (x; J, r ) + (V x t d (x; F)) T -(x r -I) + l/ 2(x r - x) T • A JtJr f d (x; ?, r ) • (x r - x) . 

Claim 9 (original): Method according to claim 8, wherein the specular 
reflection point x r (p) is determined along the length of the said reflector such that the 

diffraction travel time at the said specular reflection point x r (p) is stationary, according 
to the equation: 

p T A(V x t d (x; s , 7) + (A x J d (x; J, 7) • (x r (p) - x)) = 0 . 
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Claim 10 (previously presented): Method according to claim 8, wherein the 
specular reflection point x r and the reflection travel time t r {x r (p)\sj) are determined 
according to the following expressions: 

3c r (p) = 3f-M.F- 1 b 

t r ( x r (p); s , r ) = t d (x\ I, r ) - 1 / 2 b r - F b 
where: 

- M is a (3 x 2) matrix described by two vectors extending along the length of 
the reflector, and therefore perpendicular to the dip vector p ; 

- b is a (2 x 1) vector of first order derivatives of the diffraction travel time 
along the. reflection plane: b = M T ■ (V x t d )\ 

- F is a (2 x 2) matrix of second order derivatives of the diffraction travel time 
along the reflection plane: F= M T • (A^ x t d ) • M . 

Claim 11 (previously presented): Method according to claim 10, wherein the 
determination step uses isochronic migration maps t d {x\sj) specified for each (source 

s , receiver r ) pair involved in the migration at each image point x in the migrated 3D 
domain. 

Claim 12 (original): Method according to any one of the preceding claims, 
wherein the seismic illumination fold \{x,p) in the migrated 3D domain is estimated 
during the Kirchoff summation migration of seismic data recorded during the 3D seismic 
prospecting. 

Claim 13 (previously presented): Method for correction of seismic data 
amplitudes recorded during 3D seismic prospecting in order to compensate for the effect 
of non-uniform illumination of sub soil reflectors, comprising the steps of: 

- estimating the illumination fold I(3c, p) using the method according to claim 

1, 

- using the inverse I" 1 (x, p) of the said ratio as a weighting factor to be applied 
to each of the said seismic data amplitudes. 
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Claim 14 (previously presented): Method for selection of an acquisition 
geometry among a plurality of acquisition geometries as a function of the target of 3D 
seismic prospecting, comprising the steps of: 

- determining the illumination fold I(jc, p) by the method according to claim 1, 
for each of the acquisition geometries considered, 

- selecting the acquisition geometry providing the optimum illumination fold as 
a function of the target. 

Claim 15 (previously presented): Method according to claim 2, wherein, during 
the incrementing step, the illumination fold 1(1, p;s,r) is incremented using an 
increment function i (t d j r \s,r) according to I(3c,/?) = l(x 7 p) + i (t d J r ;s J r), the said 
increment function taking account of the difference between the diffraction travel time 
t d (3c; s, r) and the reflection travel time t r (3c r (p)\ s, r ) . 

Claim 16 (previously presented): Method according to claim 15, in which an a 
priori correction w (3c,?,r) of the illumination fold is taken into account by migration, 
comprising the step of incrementing the illumination fold I(3c, p\s,r) related to a (source 
s , receiver r ) pair by i (t d , t r ; s , r ) . w( J; s, r) . 

Claim 17 (previously presented): Method according to claim 16, wherein the 
determination step includes the second order Taylor series development of the diffraction 
travel time (x; s, r) around the image point x: 

t d (x; s, r) = t d {x\ s, r) + ( (x\ 5, r )) T • ( x r - x) + 1 / 2( x r -x) J • A xx t d (3c; s, r) • (x r - x) . 

Claim 18 (previously presented): Method according to claim 17, wherein the 
specular reflection point x r ( p ) is determined along the length of the said reflector such 
that the diffraction travel time at the said specular reflection point x r ( p ) is stationary, 
according to the equation: 

p T A(V x t d (x; s, r) + (A x J d (x; s, r) • (x r (p)-x)) = 0. 
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Claim 19 (previously presented): Method according to claim 18, wherein the 
specular reflection point x r and the reflection travel time t r {x r {p)\sj) are determined 
according to the following expressions: 

x T (p)- x -M.V~\b 

t r (x r (py,JJ) = t d (x;sJ)-\/2b r F- ] b 

where: 

- M is a (3 x 2) matrix described by two vectors extending along the length of the 
reflector, and therefore perpendicular to the dip vector p ; 

- b is a (2 x 1) vector of first order derivatives of the diffraction travel time along 
the reflection plane: b = M T • (b = M T • (V x t d ); 

- F is a (2 x 2) matrix of second order derivatives of the diffraction travel time 
along the reflection plane: F = M T • (A x x t d ) • M . 

Claim 20 (previously presented): Method according to claim 19, wherein the 
determination step uses isochronic migration maps t d (3c; J, r ) specified for each (source 

s , receiver F) pair involved in the migration at each image point x in the migrated 3D 
domain. 

Claim 21 (previously presented): Method according to claim 20, wherein the 
seismic illumination fold I(J, p) in the migrated 3D domain is estimated during the 
Kirchoff summation migration of seismic data recorded during the 3D seismic 
prospecting. 
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